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ABSTRACT 

We separate the extragalactic radio source population above ^50 ^Jy into active galac¬ 
tic nuclei (AGN) and star-forming sources. The primary method of our approach is to 
fit the infrared spectral energy distributions (SEDs), constructed using Spitzer/IKAC 
and MIPS and ITersc/iel/SPIRE photometry, of 380 radio sources in the Extended 
Chandra Deep Field-South. From the fitted SEDs, we determine the relative AGN 
and star-forming contributions to their infrared emission. With the inclusion of other 
AGN diagnostics such as X-ray luminosity, Spitzer/IKAC colours, radio spectral index 
and the ratio of star-forming total infrared flux to ^-corrected 1.4 GHz flux density, giR, 
we determine whether the radio emission in these sources is powered by star formation 
or by an AGN. The majority of these radio sources (60 per cent) show the signature 
of an AGN at some wavelength. Of the sources with AGN signatures, 58 per cent 
are hybrid systems for which the radio emission is being powered by star formation. 
This implies that radio sources which have likely been selected on their star formation 
have a high AGN fraction. Below a 1.4 GHz flux density of 1 mJy, along with finding 
a strong contribution to the source counts from pure star-forming sources, we find 
that hybrid sources constitute 20—65 per cent of the sources. This result suggests that 
hybrid sources have a significant contribution, along with sources that do not host a 
detectable AGN, to the observed flattening of the source counts at ^ImJy for the 
extragalactic radio source population. 

Key words: galaxies: active—galaxies: nuclei—galaxies: ISM—infrared: 

galaxies—radio continuum: galaxies 


1 INTRODUCTION 
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In the time since the radio sky was first imaged, it has 
been found that extragalactic radio sources can be char¬ 
acterised by two distinct types. The first is where the 
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production of the radio emission is a result of the ac¬ 
cretion of material onto a super m assive black hole at 
the centre of a galaxy (H ee J ll984ll . an active galactic 
nucleus (AGN). The second is where the emissio n stems 
from the production of stars on a massive scale llCondonI 
Il992h . a star-forming galaxy (SFG). The bright end of the 
Euclidean-normalised source counts (151MHz flux densi¬ 
ties ~10£mJyj_Jsdommated by powerful radio-loud AGN 
(e.g. IWillott et alJ l2002f) . The distinction between radio- 
loud and radio-quiet AGN was original ly based on their 
5 GH z to B band flux density ratios llKellermann et ahl 
Il989l l. Extending down to fainter flux densities, the source 
counts distribution is well described by a power law to a 
flux density of ~lmjy, whe re the distribution then flattens 
fe.g. IWindhorst et al.|[l99al . The exact cause of this feature 
is unclear. There is evidence to suggest that the flatten¬ 
ing is the result of the emergence of the SFG population, 
which therefore dominates the source counts below ~lmjy 
(lBenn_et_^ 19931: Hopkins et al.l[l9^ : ISevmour et al.l[2004 


iGondon et al. 2012li . However, there is increasing evidence 
that radio-quiet AGN could also ha ve an important contri¬ 
bution at these faint flux densities dGruppioni et al.l 199£; 


Giliegi et al.ll2003l: Simpso n et akll iooH: Seymour et al.1 200S 


Smolcic et al.11^081 : iPadovani et al.l 2009l l and therefore be 

partly responsible for the flattening. 

To investigate the nature of the extragalactic radio pop¬ 
ulation, observations at multiple wavelengths are required. 
The Extended Chandra Deep Field-South (ECDF-S) is a 
region that has been exte nsively observed a t many wave¬ 
lengths, e.g. in the radi o I Norris_et_aJj_ 200^J|^o nso et al.1 
l200d: [Kellermann et al.l l2008l: iMiller ei al.l l2008ll. infrared 


jLonsdale_et_^ 20031 : Olsen etjJJ gOOd Tavjor^t_a^ 200^ 

Dammjt^^yj201lj)^ optical (lBi 2 L.e<L^ fjgOOjjLeFJyre’et all 


20051 : Hildebrandt_et_^d2l22P^Ll22£^l^iii2ii^^kjil 20ld) and 

ii; 


X-ray dGiacconi et al.l 20021 : Lehmer et al. 20051 : Luo et ahl 
l2008i:IXue et al.ll2011^ and is therefore an ideal field for such 
a study. A vital component of a study such as this is the 
correct separation of emission from AGN and star-forming 
processes. This can be difficult because star formation di¬ 
agnostics typically either reside at wavelengths where they 
suffer from obscuration (e.g. ultraviolet and optical) and/or 
contamination from AGN emission (e.g. X-ray, mid-infrared 
and radi o). However, at far-infrared wavelengths, AGN can 
be weak llHatziminaoglou et al.l[^010lf and the thermal radi¬ 
ation from the cool dust associated with star formation can 
dominate. With the advent of the Spectral and Photometric 
Imaging REceiver (SPIRE, iGrifSn et al.1 l2010li instrument 
onbo ard the Herschel Space Observator^ ( Pilbratt et ahl 
I 2 OIOII , we now have a clear view of the Universe at these 
wavelengths, enabling us to probe the obscured star-forming 
properties of galaxies at various cosmic epochs. With the 
addition of Hersc/ief/SPIRE data to Spitzer near- and mid- 
infrared data, infrared (IR) spectral energy distributions 
(SE Ds) have been constructed for both IR-luminous AGN 
(e.g iHatzimi naoglou et al.ll2010l: ISevmour et al.ll201ll ') and 


SFGs (e.g. Brisbin_et_^ 201^ Rowa,n-Robinson et alji2010l : 

iBuat et al.ll20ld : ISvmeonidis et al.1 12 OI 3 I L Through multi- 


^ Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and 
with important participation from NASA. 


component SED fitting of AGN/SFG composite objects, 
previous studies have disentangled the contributions from 
these two proces ses to the IR emissi o n of extragalactic 
populations (e.g. iMullanev et al.1 12OIII : iFeltre et al'l I2OI3I : 
iHuang et al.ll2014 ). IR SED fitting can therefore be a pow¬ 
erful tool for detecting AGN, even in cases where they do 
not dominate the emission. 

It is now well known that there is a tight corre¬ 
lation between the far-IR and radio emission produced 
from star formation (the far-IR-Radio Gorrelation; FIRRC, 
Ide Jong et al.|[l985l : iHelou et al.1 1 19851 ). A correlation with 
such small scatter, spanning many orders of magnitude in 
luminosity, may be surprising considering that the two forms 
of emission are generated by different processes. On one 
hand the IR emission is produced by the dust in a galaxy 
that absorbs the ultraviolet light from newly-formed mas¬ 
sive stars and re-emits this energy as thermal radiation, 
while on the other, the radio emission is synchrotron ra¬ 
diation produced by the cosmic rays from supernova rem¬ 
nants. The reason for the correlation is that, as the life-cycle 
of massive stars is relatively short (~10^yrs), the two pro¬ 
cesses can be observed simultaneously i f star formation oc - 
curs continuously on longer timescales (llvison et al JH). 
Many studies have calculated the logarithmic ratio of IR 
and radio luminosit i es, qm, at both lo w redshift (2 < 0.5, 
iGondon et ah 1991 J^vis_et^ 2010h and high redshift 
(up to 2 ~3, Saiina et al.| 20081: Ibar et al.ll 2008l: llvison et al.1 


I2OI0I : IChapman et al.l 2010l : Seymour et al.ll201ll L either for 
its use in future studies or as a way of validating their re¬ 
sults. The FIRRG can be a useful tool because through radio 
observations, the IR star-forming properties of a sample can 
be inferred (or vice versa). 

In this paper, we divide the radio source population into 
sources whose radio emission is powered by an AGN or by 
star formation, partly through fitting the IR SEDs of 380 
radio sources in the ECDF-S. Herschel/SPIRE photometry, 
provided b y the Herschel Mu lti-tiered Extragalactic Survey 
IHerMES. [Oliver et al.l[20ll) . is included in the SEDs and 
we take advantage of the excellent multi-wavelength cover¬ 
age of the field. From the fitted-SEDs, we constrain the AGN 
and galaxy interstellar medium (ISM) contributions to the 
IR emission and attempt to determine if the radio emission 
is powered by star formation or by an AGN. We then cal¬ 
culate the differential radio source counts for the ECDF-S 
by source type. The paper is structured as follows. In Sec¬ 
tion (2] we describe the master catalogue, multi-wavelength 
data and the sample selection. In Section [3 we descibe the 
SED modelling and the fitting procedure. In Section [T] we 
present the results and discuss them in Section [5] Conclu¬ 
sions are presented in Sectional Throughout, we assume a 
A Cold Dark Matter Universe, with Ho = 70kms“^ Mpc“^, 
Da = 0.7 and Dm = 0.3. 


2 THE SAMPLE 

2.1 Master catalogue and photometric data 

The master source catalogue comes from the second data 
release (DR2) from the National Radio Astronomy Ob¬ 
servatory (NRAO) Very Lar ge Array (VLA) o bservations 
of the ECDF-S at 1.4 GHz llMiller et al.11201^ . The im¬ 
age has a typical sensitivity of 7.4 /ijy per 2.8 by 1.6 arcsec^ 


































































































































HerMES: Disentangling AGN and SF in radio sources 3 


beam and an area of 0.32 deg^. The catalogue contains 883 
sources above 5 a, with positional accuracies of 0.2 arcsec 
and 0.3 arcsec in right acension and declination, respectively, 
for a 5 (t source. Regarding the 1.4 GHz flux densities, 5'i.4, 
of the catalogue, it has been stated that the true radio flux 
density of an unresolved source is better represented by its 
peak flux density rather tha n by its integrated flux den¬ 
sity ||Owen_^_^orris^ 200^ . In line with other radio sur¬ 
veys, 2013) recommended the use of the peak 

flux densities for unresolved sourc es and the in t egrate d flux 
densities for resolved sources. See lMiller et al.l ll2013h for a 
description of how the radio sources were identified as re¬ 
solved or unresolved. Through the us e of multi-wavelengt h 
data and counterp art identifications llBonzini et al.l[2012ll . 
I Miller et al.l (120131) flagged 17 multiple-component sources 
likely to be powerful radio galaxies. For these sources, we 
used the summed flux densities of the individual compo¬ 
nents from the main catalogue. We also used the radio po¬ 
sitions from the main catalogue and these corresponded to 
the sources that are likely to be the core objects. For the 
multiple-component sources for which the core is not visi¬ 
ble, the average of the other component positions were given 
in the catalogue. 


Some of the data used to obtain the IR SEDs come 
from the seven-band imaging of the CDFS carried out with 
the IRAC and MIPS Spitzer instruments. To obtain IRAC 
flux densities we used the IRAC/MUSYC Public Legacy 
Survey in the ECDF-S (SIMPLE. iDamen et al.ll201lll ver¬ 
sion 3.0 catalogue. The survey has sensitivities of 1.1, 1.3, 
6.3 and 7.6 pJy, at 3.6, 4.5, 5.8, and 8.0/rm, respectively, 
at 5 ct. We cross-correlated the master radio catalogue with 
the SIMPLE catalogue using a 1.5 arcsec search radius and 
found IRAC counterparts for 770 sources (out of 883). We 
deemed the radio sources without an IRAC counterpart to 
be undetected in the IRAC bands if its position was within 
the coverage of the SIMPLE survey. See Section 13.11 for 
how we dealt with non-detections in the IRAC bands. For 
flux densities at 24 ^m, we used the catalogue produced by 
the Spitzer far-infrared Deep Extragalactic L egacy (FIDEL) 
Survey of the ECDF-S using IRAC priors dMaenelli et al.l 
l2009^ . iMagnelli et al.l (120091 1 reported the sensitivity of the 
FIDEL map as 70 /rJy at 5 cr and they presented a catalogue 
of ~9400 sources above 3 a. Taking advantage of the fact 
that the FIDEL sources have already been assigned IRAC 
counterparts, we matched the IRAC SIMPLE counterparts 
to the FIDEL catalogue using the identical FIDEL IRAC po¬ 
sitions. This gave 602 matches. For sources without 24 pm 
counterparts, we used 24 pm flux densities that were de- 
rived from the FIDEL image via the HerMES XID process 
ijRoseboom et al1l2010l l. For radio sources with IRAC coun¬ 
terparts, flux densities were extracted at the IRAC posi¬ 
tions, otherwise the radio positions were used as priors. For 
the radio sources with a 24pm counterpart, the XID 24pm 
flux densities were co nsistent with the flux densities from the 
iMagnelli et al.l ll2009ll catalogue. Furthermore, the mean flux 
density derived from the XID process of 0. 50 mjy is com¬ 
parab le to the mean flux density from the IMagnelli et al.l 
(I 2 OO 9 II catalogue of 0.51 mJy. 

The remaining photometry comes from the MIPS 70 pm 
and SPIRE 250, 350 and 500 pm bands. The 70 pm and 


SPIRE0 catalogues were generated by extracting sources 
from the 70 pm FIDEL and HerMES maps of the CDF- 
S at the 24 ^m positions via the XID process. For 
the radio sour c es wi thout a 24 ^m counterpart from the 
IMagnelli et al.l (l2009f) catalogue, the radio positions were 
used as priors. The noise in the 250, 350 and 500 pm SPIRE 
bands is dominated by confusion noise and the images 
have 1 a sensitivitie s of 5.8, 6.3 and 6.8 mJy, respectively 
liNguven et al.l [ 2 OIOII . The estimated uncertainties on the 
SPIRE flux densities take confusion noise into account. 


2.2 Sample selection and redshifts 

For consistency, we initially required coverage across all nine 
bands; IRAC 3.6, 4.5, 5.8, and 8.0 pm, MIPS 24 and 70 pm 
and SPIRE 250, 350 and 500 pm. We therefore rejected the 
radio sources from the master catalogue that fell outside 
the IRAC, 24 pm and/or SPIRE coverage; this removed 174 
sources. We made further cuts by rejecting radio sources 
that had an IRAC counterpart that was flagged in the SIM¬ 
PLE catalogue as being a Galactic star or had a flux density 
that was flagged as being contaminated by a Galactic star; 
this applied to 34 sources. We also rejected two radio sources 
that were matched to the same IRAC source and we rejected 
any radio sources that had a FIDEL 24 pm source between 
1.5 and 2.5 arcsec. This last rejection was on the basis that 
neighbouring sources are likely to be contaminating the far- 
IR flux densities. This removed 27 sources. After all these 
rejections we were left with a sample of 646 radio sources 
(hereafter called ‘full-sample’), each with photometry, in¬ 
cluding upper limits, across all nine IR bands. To estimate 
the probability of the radio sources being matched randomly 
to the counterpart sources, we shifted the radio positions of 
the master catalogue in right ascension by ± 10 arcsec, while 
keeping the declination constant and then vice-versa. We 
cross-correlated the master catalogue again with the SIM¬ 
PLE catalogue using the same search radius as used initially, 
this time with the four new positions. We found that the 
mean number of false matches for the four shifted positions 
is 83, 10 per cent of the matches for the true radio positions. 
We also did this with the 24 pm catalogue, again using the 
same search radii, and found 16 false matches (3 per cent). 

To fit model S EDs to the observed photometry, red- 
shifts were necessary. iBonzini et al. lIMi) identihed the op¬ 
tical, IR and X-ray counterparts of the radio sources in the 
VLA ECDF-S DR2 catalogue. See their Table 4 for a list of 
the optical and IR catalogues. Their X-ray data come from 
the Chandra X-ray observat ory 4 Ms survey of the central 
part of the CDFS (0.1 deg^. IXue et al.ll201ll l and a mosaic 
of four 250ks obser v ations of the extend ed field (0.28 deg^, 
iLehmer et al.lf^05l) . iBonzini et al.l ll2012f) also acquired op¬ 
tical spectra for 13 radio sources which they combined with 
data from the literature to gain redshift information for as 
many of these objects as possible. For radio sources with 
both spectroscopic and photometric redshifts, if the spec¬ 
troscopic redshift had a low quality flag, the photometric 
redshift was preferred bv iBonzini et al.1 1I2OI2I) . They found 


^ The data presented in this paper will be released 
through the Herschel Database in Marseille HeDaM 
(hedam.oamp.fr/HerMES). 
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the accuracy of the photometric redshifts to be ~6 per cent. 
All sources with a low quality spectroscopic redshift also had 
a photometric redshift and so no source had a final redshift 
that is a low-quality spectroscopic redshift. Of the 646 radio 
sources in the full-sample, 567 (88 per cent) have redshifts, 
of which 220 (39 per cent) are spectroscopic and 347 (61 per 
cent) are photometric. 

When it came to htting the IR SEDs of the radio 
sources, a rest-wavelength range of 2—1000 pm was used. 
For increasingly high—redshift sources, an increasing num¬ 
ber of the IRAC bands are redshifted outside this wavelength 
range. This decreases the number of photometric points used 
during the htting. To ensure there was a sufEcient number 
of degrees of freedom to constrain the AGN and galaxy ISM 
SED components, we therefore also rejected sources above 
z ~ 1.4. This removed 123 sources. Next the IR SEDs of the 
remaining 444 sources were htted and the effect of redshift 
uncertainties for the sources with photometric redshifts was 
investigated (see Section r3.1l for a full description of the SED 
htting procedure). We htted the SED of each source three 
times; first using the redshift measurement and then using 
the redshift 3 a Gaussian upper and lower bounds. Each 
best-ht returned a minimum value and the difference in 
these values from the best case to the worse case scenar¬ 
ios was measured. If the difference in x^ was > 9, the 3 a 
conhdence level for one ‘interesting par ameter’ when usin g 
a standard x? minimisation technique llPress et al. 
then the hts were not consistent at a signihcance of 3 (t. This 
meant that the redshift uncertainty was too large to con¬ 
strain the AGN and galaxy ISM SED components for any 
source that met this Ax^ criterion. Here, the term ISM refers 
to the IR emission from a galaxy that is attributed to star 
formation i.e. cold dust and polycyclic aromatic hydrocar¬ 
bon (PAH) emission. This criterion occured for 64 sources. 
We were therefore able to separate the AGN and galaxy ISM 
contributions to the IR emission for 380 of the 646 sources 
in the full-sample. These 380 radio sources constitute a sub¬ 
group of the full-sample (hereafter called ‘sub-sample’). The 
majority of the work presented here is based on this sub¬ 
sample, although we do consider the full-sample when cal¬ 
culating the radio source counts. The redshift distribution of 
the sub-sample peaks at a ~0.6 and is displayed in Figure [T] 
There are non-detections for four sources in the four IRAG 
bands and 26 sources at 24 pm. There are X-ray counter¬ 
parts for 130 sources. The full- and sub-sample selections 
are summarised in Table [T] 


2.3 Ancillary data 

Radio spectral indices were desired for the sub-sample and so 
we also determined 5.5 GHz hux densities. To do this we used 
the catalogue of 145 sources from th e ATLAS 5.5 GHz survey 
of the EDGF-S dHuvnh et al.l[201^ L The ATLAS image has 
an almost uniform sensitivity of ~12pJy rms, a beam size 
of 4.9arcsec x 2 arcsec and an area of 0.25 deg^. We cross- 
matched the master catalogue to the ATLAS catalogue us¬ 
ing the radio positions and a 1.5 arcsec search radius. There 
were ATLAS counterparts for 60 sources in the sub-sample. 
As with the master catalogue, for unresolved sources we used 
the peak hux densities and corresponding uncertainties for 
the 5.5 GHz hux densities. For resolved sources, we used the 
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Figure 1. The redshift distribution of the sources in the sub¬ 
sample (for definition see text). The black solid, red dashed and 
blue dot-dashed lines show the total, spectroscopic and photo¬ 
metric redshift distributions, respectively. Also shown as a green 
dotted line is the redshift distribution of the sources excluded 
from the sub-sample because their photometric redshift uncer¬ 
tainties were too large to constrain the AGN and galaxy ISM 
SED components, see Section |3. II 


Table 1. Breakdown of sources in the master radio catalogue 
(see Section ESI. The full-sample refers to the sources from the 
master radio catalogue for which there was photometry, including 
upper limits, across all nine IR bands used in this study. The sub¬ 
sample refers to the sources in the full-sample for which we were 
able to separate the AGN and galaxy ISM contributions to the 
IR emission (see Section [3T}. 


Number of sources in master catalogue 833 

Number of sources in full—sample 646 

Number of sources in full—sample with redshifts 567 

Number of sources in sub—sample 380 

Number of sources in sub—sample with X-ray data 130 


integrated hux densities and assigned fractional uncertain¬ 
ties equal to the peak hux density fractional uncertainties. 

There were ten multiple-component sources in the AT¬ 
LAS catalogue and the individual components were listed 
separately. As only a single central position was given for the 
multiple-component VLA sources, there were larger ohsets 
than expected between the extended VLA sources and their 
ATLAS counterparts. Except for two cases, all multiple- 
component ATLAS sources were matched to multiple- 
component VLA sources. For one of these cases, the AT¬ 
LAS source was interpreted to be two separate sources in 
the VLA catalogue and for the other case, because of the 
extended nature of the ATLAS source components, there 
was a separation larger than 1.5 arcsec between the AT¬ 
LAS component positions and the VLA position. For the 
other eight multiple-component ATLAS sources, the individ¬ 
ual ATLAS component hux densities were combined to give 
total 5.5 GHz hux densities. The radio spectral index, ctrad, 
dehned by 5"^ oc , was then calculated for the 60 sources 
with 5.5 GHz counterparts using 1.4 and 5.5 GHz hux den- 
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sities. For the sources in the sub-sample without a 5.5 GHz 
counterpart, we assumed a radio spectral index of —0.8. 

F or th e sources with X-ray data from the IXue et aP 
ll201lll and iLehmer et al.l (120051) catalogues, we calculated 
rest-frame 0.5—8keV luminosities, Lq.s-s, using observed- 
frame 0.5—8keV fluxes, /o.s-s- The equation we used was 
To. 5-8 = dvr Dl /o.5-8 (1 + where F is the X-ray pho - 

ton index and a value of 1.8 was assumed dXue et al.l 201lll . 
Note these luminosities are not absorption-corrected, but be¬ 
cause such a correction would only increase the luminosity, 
any sources identified as AGN through their X-ray luminosi¬ 
ties (see Section r3.2ll would remain identified as AGN if their 
luminosities were corrected. This could, however, mean that 
heavily-obscured AGN may not be identified through their 
X-ray emission. 


3 DETERMINING THE AGN AND ISM 
CONTRIBUTIONS TO THE RADIO 
OUTPUT 

3.1 IR SED fitting 

For the sub-sample of 380 sources, we determined the AGN 
and galaxy ISM contributions to their IR emission. To do 
this we used an IDL fitting procedure based on DecompII0 
dMullanev et al.ll201ll ~). We made certain asssumptions on 
the shapes of AGN and ISM SEDs though our choice of tem- 
pla tes. We cha r acteri sed the emission from the AGN using 
the lElvis et al.l dl994l i radio-quiet QSO SED and the Spitzer 
Wide- Area Infrared Extr agalactic (SWIRE) template li¬ 
brary dPolletta et al.l[2007l l . We used five of the seven empir¬ 
ical AGN templates from the library and excluded the two 
Seyfert galaxy SEDs. These two Seyfert SEDs exhibit strong 
PAH features and therefore there is a non-negligible contri¬ 
bution from star formation in the h ost galaxies. To model the 
emissi on from the ISM we used the lSiebenmorgen fc Krfigell 
d2007l. hereafter SK07) library of SEDs for starburst and 
ultraluminous IR galaxies. These were created using radia¬ 
tive transfer models and the library contains ~7000 model 
SEDs that vary in total luminosity, starburst nucleus ra¬ 
dius, visual extinction, ratio of the luminosity of OB stars 
embedded in dense clo uds to the total luminosi ty and dust 
density in these clouds. ISvmeonidis et akl (l2013l l found that 
up to 2 ~2, some of the SK07 SEDs were not representa¬ 
tive of IR-luminous (log(I/iR/LQ) > 10) SFGs. They char¬ 
acterised the SK07 SEDs in terms of a ‘flux’ parameter 
which was defined as Ltot where Ttot is the lumi¬ 

nosity of the SED in units of Lq and R is the nuclear star- 
burst radius in units of kpc. They found that the SEDs with 
log(fiux / Lq kpc“^) <8 or log(flux / Lq kpc“^) >11 were 
not good fits to their sample. Taking their results into ac¬ 
count, we therefore considered only the SK07 SEDs with a 
flux in the range 8 < log(flux / Lq kpc“^) <11. Of the ~7000 
SEDs in the original library, 4765 met this criterion. By only 
using these SEDs as the templates for our ISM component, 
we have assumed that the IR ISM SEDs of our sub-sample 
are similar to those of IR-luminous SFGs. If the IR emission 
of our sources is star-formation-dominated then this would 
imply they are also IR-luminous SFGs and this assumption 

® https://sites.google.com/site/decompir 


is therefore reasonable. If, on the other hand, these sources 
are quiescently forming stars then they could have ISM 
SEDs that are different to those of IR-luminous SFGs. They 
could perhaps ha ve SEDs more like the cirrus -dominated 
SEDs observed by iRowan- Robinson et akl (l201Cll l. However, 
if they are quiescently forming stars then the AGN will likely 
be the dominant component in their near- to mid-IR SED. 
This then means that the disputed nature of an ISM com¬ 
ponent that has only a small contribution would not be a 
major issue. For consistency, we therefore fitted the 4765 
SK07 SEDs to all sources in our sub-sample. 

At near-IR wavelengths, there is a contribution to the 
total output of a galaxy from direct starlight. Although di¬ 
rect starlight is taken into account in the radiative trans¬ 
fer models for the SK07 templates, we desired an indepen¬ 
dent starlight component to increase the flexibility of our 
SED-fitting model. We therefore included a third compo¬ 
nent to repres ent this emission. For th is component, we 
made use of the lBruzual fc Gharlotl (120031 ') stellar population 
(SP) library. The library contains 78 templates which vary 
in star-forming history, stellar age (from 5 Myr to 12 Gyr) 
and metallicity (from Z = 0.008 to 0.05) and cover the wave¬ 
length range from 9 nm to 160 /rm. While the SEDs vary 
greatly at the shorter wavelengths, longward of ~2 fim they 
are ah well approximated by a Rayligh-Jeans tail. There¬ 
fore, by only htting at wavelengths >2/im (rest-frame), we 
were able to use just one SP template as the starlight com¬ 
ponent. With the addition of this third component to our 
model, the direct light from the SK07 templates needed to 
be removed. This was done by subtracting the SP template, 
having hrst normalised it to the SK07 template flux at 2 fim. 
These modihed SK07 templates then purely represent the 
emission from the ISM of the galaxy. The normalisation of 
the SP template was hxed only when removing the direct 
starlight from each SK07 template. The SP component nor¬ 
malisation was a free parameter during the SED htting. 

To account for the reddening of the AGN that can occur 
in some sources we made an assumption that the dusty envi¬ 
ronments of high-reds hift AGN can be modelled by a Galac¬ 
tic extinction curve (Saiina et al.l 20071: iLacv et al.1 20071: 
Sevrnour et al.1 l2008l : Drouart_et_^ 20121 : Rawlings et al.1 


2 OI 3 II . We chose the |d rain d (|2Q0j) extinction curve that 


uses Rv =3.1, the typical value for the Milky Way Galaxy, 
where Rv is the ratio of total extinction, Av, to selective ex¬ 
tinction, Ab — Av , in the U band (i.e. Rv ~ Ay/ (Ab — Av)) 
and we applied this curve to the AGN component only. Us¬ 
ing such a curve could result in an AGN component with a 
large silicate absorption feature. This approach assumes that 
the obscuring medium acts as a simple foreground screen 
and is close to the AGN but not located i n the torus it¬ 
self. T hrough their clumpy torus modelling, iNenkova et al.1 
ll 2008 tl showed that the silicate absorption strength would be 
weak (r 9 , 7 (/m < lj3 T dust causing the AGN reddening 
was situated in the torus. Therefore, to also consider this sce¬ 
nario, we reapplied the extinction curve with the silicate fea¬ 
ture removed by interpolating over the relevant wavelength 
region. This allowed for the case of an AGN component with 


^ Throughout the paper 79,7 nm represents the 9.7/rm apparent 
optical depth and is defined as the natural log of the ratio of 
continuum emission to observed emission. 
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a large amount of reddening but no silicate feature. There 
was no need to apply extinction to the ISM component be¬ 
cause extinction was already applied to the models used to 
create the SK07 templates. We did not apply extinction to 
the direct starlight component because this would decrease 
the slope of the template and this component would start 
to become degenerate with the AGN component over the 
wavelength region in which the starlight component most in¬ 
fluences the SED-fitting (A < 10/rm rest-frame). This would 
have caused the direct starlight component to become re¬ 
dundant. 

During the IR SED fitting we used six free parameters: 
(i) the chosen ISM template, (ii) the ISM component 
normalisation, (iii) the chosen AGN template, (iv) the AGN 
component normalisation, (v) the direct starlight compo¬ 
nent normalisation and (vi) the extinction of the AGN 
component i.e. A(A). An upper limit restriction was placed 
on A(A) such that Av ~ 30 mag. This limit corresponds to 
a large 9.7 pm silicate absorption depth of for 

the extinction curve used. A silicate absorption depth of 
this size has been observed in galaxies selected to be dust y 
i.e. large mid-IR-optical colours (e.g. ISaiina et al.l l2007ll . 
Only six free parameters have been listed here and so we 
make the caveat that there are other parameters embedded 
within the ‘choice of ISM template’ parameter. These are 
the free parameters of the radiative transfer modelling 
which produced the SK07 SED library. We fitted all three 
components of our model in our rest-wavelength range and 
the expression for the model, Sv, in mjy was 

S'.(A) = fcatt^A) ^ + k,cRX) (1) 

where fca and av{X) are the normalisation factor and flux 
density of the AGN component respectively, fcb and b^{X) 
are the normalisation factor and flux density of the ISM 
component and kc and Ci/(A) are the normalisation factor 
and flux density of the stellar population component. For 
each ISM template, minimisation was performed on the 
SED, simultaneously solving for fca, fcb, fee and A(A). The 
photometric data were assigned a minimum uncertainty of 
10 per cent so that individual data points were not weighted 
too high in the estimation during the fitting. The mini¬ 
mum error is justified by the need to take into account the 
limitations of the templates which are being used in the 
modelling. Both the AGN and the starburst model compo¬ 
nents are taken from libraries of discrete SEDs rather than 
continuous parameter-based models, and the only degrees 
of freedom allowed beyond the choice of the library tem¬ 
plate are the normalisations of the templates and the degree 
of dust reddening of the AGN. SK07 compare their library 
SEDs to a number of well know star-forming galaxies, and 
show the residuals of their fits in their figure 7. The scatter 
shown in the various panels of that plot typically exceeds 10 
per cent, though the statistical and systematic uncertain¬ 
ties on the data points contribute as well as the limitations 
of the models. Bearing that in mind, 10 per cent appears 
to be a reasonable estimate for the data-model mismatch 
which is likely to come from the limited set of SK07 SED 
shapes. In practical terms, such an approach is required be¬ 
cause otherwise the statistical weight of the near-IR bands 
would be much greater than the far-IR bands, and the fitting 


of the near-IR photometry would dominate the goodness 
of fit. Without the adjustment to the photometric uncer¬ 
tainties at short wavelengths, we found that the best-fitting 
templates often fail to reproduce the far-IR emission, and 
therefore the SED fitting process would not provide reli¬ 
able constraints on the ISM component. We note that in 
adopting a minimum error of 10 per cent we are fol l owing 
the same approach as t hat taken bv ISvmeonidis et al.l (I200S ) 
and lSvmeonidij ll2008l l to model the SEDs of infrared galax¬ 
ies using the SK07 templates. 

We ensured that source non-detections in the IRAC 
bands had less weight than detections during the SED fit¬ 
ting. Our technique was to assign flux densities and 1 a un¬ 
certainties that were equal to half the 5 a survey limits of 
the SIMPLE survey, giving a S/N = 1. For example, sources 
not detected at 8 /rm were given an 8 fim flux density and 
uncertainty of 7.6 /iJy/2 = 3.8 /rJy. With this technique, at 
each wavelength in the 3.6—8 /rm range, non-detections had 
individual x^ values of < 1 for fitted models that predicted 
a flux density below the survey sensitivity. This meant such 
a model was acceptable at the 1 a level at each wavelength 
where there was a non-detection. If a model predicted a flux 
density above the survey sensitivity then a non-detection 
would have an individual x^ value of > 1 and the model 
would not be acceptable at the 1 a level. 

To determine the confidence in the best-fitting ISM 
component normalisation (feb), we found the extent to 
which it can vary until the condition x^ ~ Xmin = 
is met, letting a = 1 or 3 to find the 1 and 3 cr x^ confidence 
levels for one ‘interesting parameter’. During this process 
the free parameters: fea, fee and A(A) were allowed to vary 
also. Total IR (8—lOOO/rm) star-forming fluxes, S'ir,sf, 
were derived from the best-fitting ISM components. These 
were co nverted into lumi nosities, Tir.sf and then into SFRs 
via the Ik ennicu t3 (liggsl i relation 

SFR _ Lir,sf 

IMoyr-i ^ 5.8 x lO^ Lq ^ ’ 

The r elation assumes solar abundances and a ISalpeteJ 
(Il955ll initial mass function (IMF). Uncertainties in I/ir,sf 
and the SFRs were derived from the uncertainties in the ISM 
normalisations. As A(A) was allowed to vary throughout, the 
uncertainties in the above parameters have also taken into 
account the uncertainty in the level of extinction applied. 
For example, where the extinction is not well constrained, 
the confidence levels on the above parameters include AGN 
components both with and without extinction. 

For each source, we also sought to determine if AGN 
and ISM components were required in the fitted SEDs. To 
do this we applied a more conservative significance cut than 
used previously as this determination will have important 
consequences later on for this study. Firstly, the reduction 
in the minimum x^ as a result of the inclusion of an AGN 
component with respect to a purely star-forming template 
was measured. If this Ax^ was > 25 then at the 5 a con¬ 
fidence level for one ‘interesting par ameter’ when us ing a 
standard x^ minimisation technique (IPress et al.lll992h . the 
fits with and without the component were not consistent 
at a significance of 5cr. For any source that met this Ax^ 
criterion we therefore determined that at the 5 a level, an 
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AGN component was deemed to be required in the IR SED, 
i.e. an AGN had been detected. Similarly, if the inclusion of 
an ISM component reduced the minimum by > 25 with 
respect to a model consisting of only an AGN component, 
then an ISM component was deemed to be required at 5 a. 

3.2 AGN indicators 

For the 380 radio sources in the sub-sample we aimed to 
determine two characteristics: (1) whether or not an AGN 
is present in the galaxy; and (2) whether or not the ob¬ 
served radio emission is powered by an AGN. To determine 
characteristic (1), we made use of multi-wavelength data, 
as different regions of the AGN structure can be probed by 
their emission at different wavebands. Along with fitting the 
IR SEDs with AGN and ISM components, we also applied 
four observational-based indicators to each source and iden¬ 
tified that an AGN is present if any of the following criteria 
were met: 


Next we determined characteristic (2), i.e. whether or 
not the observed radio emission is powered by an AGN, for 
the 380 sources in the sub-sample. If we have been unable to 
show that an AGN is present in the source then we inferred 
the radio emission of the source is powered by star forma¬ 
tion. If there is an indication of jet activity in the source (i.e. 
criteria iii or iv) then we inferred the radio emission of the 
source is powered by an AGN, regardless of the results from 
the other indicators. If there is an indication of black hole 
accretion (i.e. criteria i or ii or through the SED fitting) but 
there is no sign of jet activity then the source is referred to 
as a hybrid. We use this nomenclature for the latter source 
type because the use of multi-wavelength data has revealed 
that an AGN is present in the source. However, as there is 
no sign of jet activity, it is not clear whether or not the radio 
emission is powered by the AGN. 


4 RESULTS 


(i) an X-ray luminosity greater th an 10"^^ ergs~^ . _ 

(ii) IRAC colours that satisfy the iDonlev et all ll2012l l 
criteria for having an AGN-dominated near- to mid-IR SED. 

(iii) a radio spectral index with a 3 ct upper bound < — 1 
or a 3 (t lower bound > —0.5. 

(iv) if the sources resides more than 3 (t above the 
FIRRC. 


For criterion (i), an X-ray luminosity cut of 10“*^ ergs“^ 
is often used by extragala ctic studies to d i fferentiate between 
AGN and SFGs (e.g. seelFabbianolll989l : iMoran et al.lll99^ : 
iBrandt fc Hasingeill2005l ). However, this is not a theoretical 
upper limit to star-forming X-ray luminosities and so we 
make the caveat that it is possible for a small number of 
extreme star-forming sources in the sub-sample to have an 
X-ray luminosity above this cut. 

For criterion (iii), this spectral index range was used 
since it has been shown that synchrotron radiation at ra¬ 
dio wavelengths, as a result of star-forming processes, has 
a spectrum with an index in the range —1.0 to —0.5 
llThompson et all l2006l l. Althongh AGN could have radio 
spectral indices within this range, if a source has an index 
outside of this range then we infer that this is because of the 
presence of an AGN. 

For criterion (iv), we calculated gm for each source 
us ing a slightly modi fied form of the equation first defined 
by H.gIou_et_aL llOSSi) . Follo wing other studies dlvison et al.l 
l20l3 i' Chapman et al.llioi^ '). the equation we used was 


qm = logic 


>S'ir,sf/3.75 X 10^^ W m ^ 

Sl.4GHz/Wm-2Hz-l 


(3) 


Radio 1.4 GHz flux densities were fc-corrected using the ra¬ 
dio spectral indices we have derived (see Section 12.311 for 
sources with 5.5 GHz c ounterparts, or oth erwise using the 
assumed index of —0.8 (llvison et al.ll201(]| l. If the gm value 
wa s less than th e 3 a lo wer bound of the median qm derived 
by llvison et al.l ll2010ll . gia = 2.40 ± 0.24, then we inferred 
that the radio emission is produced by an AGN. If none of 
the AGN criteria (i—iv) was met and an AGN component 
was not required in the fitted SED then we have been unable 
to show that an AGN is present in the source. 


Of 380 sources in the ECDF-S with fitted IR SEDs, the ma¬ 
jority (227; 60 per cent) are identified as an AGN through 
their SED and/or by at least one of the observational-based 
indicators. For 253 of the 380 sources in the sub-sample, an 
ISM component is required in their fitted SEDs at the 5 a 
level (Ax^ >25). All but five of these sources are detected 
in at least one SPIRE band at 3 (t. For the remaining 127 
sources, an ISM component is not required at the 5 cr level 
and thus we do not significantly detect ongoing star forma¬ 
tion in these sources. For 123 of the 380 sources, an AGN 
component is required at the 5 a level. All of these sources 
are detected in at least two IRAC bands at 3 (t but only 
15 have an AGN fraction of the total IR luminosity that 
is > 50 per cent. The fitted SEDs of the 123 sources have 
a median reduced of 6-4. The ISM templates cannot re¬ 
produce the high near- to far-IR ratios observed for these 
sources and so an additional AGN component significantly 
improves the fit. For 92 sources, neither an ISM component 
nor an AGN component is required. This is a result of the 
quality of the photometric data being too low to be able to 
contrain the contribution from each component. In Figures 
[2] and 131 we show several examples of best-fitting models 
to the IR SEDs for the sub-sample. In the right-hand pan¬ 
els of Figures [2] and [S] we show examples of SEDs that do 
not require an ISM and AGN component, respectively, at 
5a. However, this does not mean these components are not 
present in the best-fitting models. 

The number of sources identified as an AGN decreases 
with increasing number of indicators considered and only 
two sources are identified as an AGN by all four indica¬ 
tors and through their SED. The indicators that find the 
highest and lowest AGN fractions are the X-ray luminos¬ 
ity (67/130; 52 per cent) and IRAG colours (15/380; 4 per 
cent) indicators, respectively. Of the 123 sources identified 
as AGN through their SED, 84 sources were not identified 
as an AGN by any of the indicators. All SEDs shown in the 
left panels of Figure [3] are examples of the fitted SEDs of 
sources for which an AGN component was required at 5 cr 
but none of the indicators identified an AGN. The results for 
the various indicators are displayed in Table [31 The results 
are displayed in matrix form which allows the results for 
different combinations of indicators to be discerned. Over- 
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Figure 2. Examples of best-fitting models, denoted by a solid 
red line, for IR SEDs that, at ^5 a, require an ISM component 
(left panels) and do not require an ISM component (right panels). 
Included are the separate components: AGN-green dashed line; 
ISM-blue dotted line; and direct starlight-orange dot-dashed line. 
Upper limits are Scr. 


Figure 3. Examples of best-fitting models, denoted by a solid red 
line, for IR SEDs that, at ^ 5 cr, require an AGN component (left 
panels) and do not require an AGN component (right panels). 
Included are the separate components: AGN-green dashed line; 
ISM-blue dotted line; and direct starlight-orange dot-dashed line. 
Upper limits are 3 a. 


Table 2. Number of sources identified as AGN by the various 
indicators. Column labels denote the indicator labels described 
in Section 13.21 


Indicator 

X-ray 

(i) 

IRAG 

(ii) 

a 

(iii) 

FIRRG 

(iv) 

SED-fitting 

X-ray 

67 

— 

— 

— 

— 

IRAC 

11 

15 

- 

- 

- 

Q 

5 

2 

13 

— 

— 

FIRRC 

23 

6 

12 

94 

- 

SED-fitting 

28 

13 

5 

18 

123 

Total 

67 

15 

13 

94 

123 


all, we find that for 153 of the 380 sources (40 per cent), 
the radio emission is powered by star formation (hereafter 
called ‘SF-powered’ sources). For 95 sources (25 per cent), 
the radio emission is powered by an AGN (hereafter called 
‘AGN-powered’ sources). The remaining 132 sources (35 per 
cent) have been found to be hybrid^ 

® See Section 13.21 for the SF-powered, AGN-powered and hybrid 
criteria. 


We nex t consider the flu xes, in the sense em¬ 
ployed by ISvmeonidis et al.) (|2013l ) to characterise 
the SK07 templates, of the best-fitting ISM com¬ 
ponents for the different types of sources. We find 
that the ISM components of the SF-powered sources 
have comparable fluxes ((log(flux/L q kpc“^)) = 8.83, 
standard error = 0.06) to those of the AGN-powered 
((log(flux/Z /0 kpc“^)) = 8.71, standard error = 0.09) and 
hybrid ((log(flux / Lq kpc“^)) = 8.81, standard error = 0.07) 
sources. Although the 380 radio sources in our sub-sample 
have been distinguished in terms of whether an AGN is 
detected (at any wavelength), we find that there is no 
apparent difference in the shapes of their IR ISM SEDs. 
Overall, the majority of sources have star-forming total IR 
luminosities that would qualify them as either LIRGs or 
normal IR galaxies (NIRGs; log(LiR/L 0 ) < 11). ULIRGs are 
less common and there are no HyLIRGs. Table [3] shows the 
breakdown of sources into the luminosity categories. Trans¬ 
lating star-forming IR luminosity into SFR, there are some 
high-redshift sources with extreme SFRs (> 1000 M 0 yr“^). 
Comparable SFRs have been o bserved in high-r e dshift 
Herschel-selected galaxies (e.g. iRodighiero et all 1201 ll : 
lElbaz et al.l I 2 OII ) and other non-Her schel-selected galax¬ 
ies (e.g. SFGs, Seymour et al.l l2008l : high-redshift radio 
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Figure 4. IRAC colours for the sub-sample. The area enclosed by the dashed lines is the region in colour space in which sources that 
have an AGN-powered near- to mid-IR SED reside, according to the lDonlev et all ll2012t) criteria. SF-powered, AGN-powered and hybrid 
sources are shown as black stars, red circles and orange diamonds, respectively. Sources with upper limits to their IRAC flux densities 
(see Section |2.Ill lie off the plot. 


Table 3. Number of NIRGs, LIRGs and ULIRGs in the sub¬ 
sample according to their star-forming total IR luminosity. 


IR luminosity range 

Type 

# 

Fraction 

[%] 

log(TlR,SF/L0) < 11 

NIRG 

162 

43 

11 ^ log(LiR,sF/L©) < 12 

LIRG 

203 

53 

12 ^ log(LiR_sF/L©) < 13 

ULIRG 

15 

4 


galaxies, Rawlings et al. I I 2013 I: and sub-millimeter galaxies, 
iPope et alfeood i. For all 380 sources, we find a median SFR 
of 20Moyr~^, with a median absolute deviation (MAD) of 
20 Mq yr“^. For the SF-powered and AGN-powered sources, 
we derive median SFRs of 40 (MAD = 25) Mq yr“^ and 
5 (MAD = 5) M© yr“^, respectively, while for the hybrids, 
we find a median SFR of 20 (MAD = 15) M© yr“^. At low 
redshift, the fraction of sources powered by star formation 
that host an AGN (i.e. hybrid/(hybrid-|-SF-powered)) is 
high: 0.7 at 2 : < 0.5. This fraction decreases at high redshifts 
to 0.4 for 0.5 < 2 < 1 and 0.3 for 1 < 2 < 1.4. 

The IRAC colour distribution of the sub-sample is 
shown in Figure [4] The AGN-powered and hybrid sources 


have a wide distribution in IRAC colo ur space and can lie 
well outside the region governed by the lDonlev et al] (I2OI2II 
criteria. In Figure O we plot rest-frame 1.4 GHz luminos¬ 
ity versus star-forming total IR luminosity to show which 
sources follow the FIRRC. The AGN-powered sources, by 
definition, lie above the FIRRC because of the additional 
AGN component to radio emission. Hybrids on the other 
hand reside close to the FIRRC in a similar manner to the 
SF-powered sources. In Figure|6]we show the gm distribution 
as a function of 1.4 GHz luminosity. The additional compo¬ 
nent to the radio luminosity for the AGN-powered sources 
results in lower gm values. In a similar manner to Table [S] 
we also signify how the AGN indicator results vary for each 
source. For instance, it is visible that many of the AGN 
identified through their X-ray emission lie on the FIRRC. 
For these particular sources we next plot their excess ra¬ 
dio luminosity above that predicted by the FIRRC as a 
function of X-ray luminosity in Figure [7] We determine the 
Spearman rank coefficient, p, and the significance of the de¬ 
viation of p from zero, cr, for the AGN-powered and hy¬ 
brid sources. We find a moderate correlation for the AGN- 
powered sources {p — 0.66, cr = 10“®) and no correlation for 
the hybrids (p = 0.07, cr = 0.55). 
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Figure 5. Rest-frame 1.4 GHz luminosity as a function of star-forming IR luminosity (lower abscissa) and SFR (upper abscissa) for 
SF-powered (top panel, black stars), AGN-powered (middle panel, red circles) and hybrid (bottom panel, orange diamonds) sources. For 
sources with a < 1.5 cr detected star-forming IR luminosity, 3 (t upper limits are shown. The solid line represents the FIRRG with a qir 
value equal to the mean SF-powered giR of 2.30, while the dashed lines represent the di 1 cr limits, cr^^sF = 0.23, of the correlation. 
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Figure 6. Value of qiji as a function of 1.4 GHz luminosity for SF-powered and hybrid sources (top panel) and AGN-powered sources 
(bottom panel). For each source, we show which indicators flagged an AGN with different coloured circles. The dashed line in the top 
panel represents the mean qiji of the SF-powered sources and the dotted line in both panels identifies the region below which sources 
are more than Scr above the FIRRC. 
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Figure 7. Excess radio luminosity above that predicted by the FIRRC as a function of X-ray luminosity for AGN-powered (red circles) 
and hybrid (orange diamonds) sources. 


Although it can be seen in Figures [5] and [6] that the 
hybrid sources have similar giR values to the SF-powered 
sources, we next investigate their gm distribution more 
quantitatively. Firstly, for four 1.4 GHz flux density bins of 
equal size in log space, the number of SF-powered sources 
per giR, bin width of 0.2 as a function of gin, is shown in 
Figure [S] The qiR, distribution of the SF-powered sources 
is close to a Gaussian distribution and these sources have 
a mean qm of 2.30 and a standard deviation of 0.23. This 
value for qm is consi stent with those fou nd by other stud- 
ies (qiR = 2A0 ± O .lO. ISevmour et al.ll201ll: qiR = 2.40 T 0.24, 
llvison et al ]|2010l : OTR = 2.54 ±0.21. ISaiina et al.ll2008l) . The 
hybrid sources have a comparable mean qiR of 2.26 and a 
standard deviation of 0.25. For five 1.4 GHz flux density bins 
of equal size in log space, the number of hybrids per qiR bin 
width of 0.2 as a function of qiR is given in Figure [9] If the 
same process is powering the radio emission in the hybrids, 
i.e. star formation, then it can be expected that their qiR 
distribution would be Gaussian also and it would be spread 
uniformly around the mean qiR of the SF-powered sources. 
On the other hand, if the radio emission is powered by an 
AGN for a significant fraction of these sources, then there 
should be a significant excess in the number of sources with 
a qiR below the mean. For each flux density bin, to measure 
any potential excess as a function of qiR, firstly a Gaussian 
function with a mean and a taken from the SF-powered 
sources in the same flux density bin is normalised to the 
observed distribution of the hybrids. The normalisation is 
such that the integral of the Gaussian curve from the mean 


SF-powered qiR out to infinity, is equal to the number of hy¬ 
brids with a qiR value above the mean SF-powered qiR. The 
normalised Gaussian curve is then integrated in each qiR 
bin range and subtracted from the binned qiR distribution 
to create a residual distribution. The right-hand panels of 
Figure [9] show that this residual distribution is small across 
all flux density bins. By summing the residual distribution 
across all qiR and flux density bins, we find that there is 
an excess of 14 more hybrids with a qiR below the mean 
SF-powered qiR than above. However, this number is not 
significant as it is comparable to the Poisson uncertainty in 
the total number of hybrids in the sub-sample (■\/132 ~ 11). 

We test the sensitivity of our results with respect to 
the FIRRG indicator by increasing the threshold such that a 
source is identified as an AGN only if it resides more than 5 a 
away from the FIRRG. The overall number of sources found 
to contain an AGN decreases from 227 to 200 sources. The 
number of SF-powered and AGN-powered sources as a result 
of this change is 180 and 53, respectively. There are also now 
147 hybrids. Across all radio flux density bins, there are now 
21 more hybrids with a qiR above the mean SF-powered qiR 
than below. With a Poisson uncertainty in the number of 
hybrids of \/147 ~ 12, this number is still not significant. 
Finally, the fractions of sources powered by star formation 
that host an AGN across the three redshift bins (z < 0.5, 
0.5 < z < 1 and 1 < z< 1.4) decrease by < 2 per cent. 

Next the Euclidean-normalised radio source counts are 
calculated for the ECDF-S, using the same radio flux density 
bins as for the qiR distributions, with the extension to higher 
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Figure 8 . The giR distributions of SF-powered sources for different 1.4 GHz flux density bins. The flux density bin ranges are as 
follows: (1) 0.030—0.067 mjy; (2) 0.067—0.150 mjy; (3) 0.150—0.333 mjy; and (4) 0.333—0.742 mjy. The dashed vertical line in each 
panel represents the mean gur value. 


flux densities to incorporate the AGN-powered sources. For 
this calculation, we include the 266 sources in the full-sample 
that have unknown redshifts or have redshift uncertainties 
that were too large to be able to constrain the AGN and 
galaxy ISM contributions (see Section 12.21 hereafter ‘no-z 
sources’). We show the Euclidean-normalised source counts 
for the ECDF-S in Figure [10] along with the counts for other 
fields as a comparison. In Figure [TT] and Table UJ we break¬ 
down the source counts into the SF-powered, AGN-powered, 
hybrid and ‘no-z’ components. Additionally in Figure 1111 
we show the LERG component of the AGN-powered source 


counts along with the models f rom the simulations for the 
Square Kilometre Array (SKA. Iwilman et akilioosh . Frac¬ 
tional uncertainties in the counts have been calculated by 
assuming a Poisson distribution and are equal to 1 / v/iV, 
where N is the number of sources in each bin. We have cor¬ 
rected for the reduced coverage of the EGDF-S observation 
at the fainter flux levels using Figure 2 from [Miller et al.l 
(l2013ll . which shows the fractional area covered by the sur¬ 
vey as a function of sensitiv ity. When creating the sensitiv¬ 
ity map, [Miller et al.l ll2013ll removed sources with peak flux 
densities greater than 150 /rjy. It is assummed that along 
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Figure 9. Left panels: The qir. distributions of hybrid sources as solid orange lines and the normalised SF-powered Gaussian distributions 
as black dashed lines. Right panels: The residual distributions as solid orange lines, after the normalised SF-powered Gaussian distributions 
have been subtracted from the binned hybrid distributions. The various panels vertically show the qi^ and residual distributions 
for different 1.4 GHz flux density bins. The flux density bin ranges are as follows: (1) 0.030—0.067 mjy; (2) 0.067—0.150 mjy; (3) 
0.150—0.333 mjy; (4) 0.333—0.742 mJy; and (5) 0.742—1.653 mJy. The dashed vertical line in each panel represents the mean qir value 
of the SF-powered sources. 
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Figure 10. The Euclidean-normalised 1.4 G Hz source counts for the radio sources i n the ECDF-S (blue squares). Also given are the 
sourc e counts for the following fields: SXFD jSimpso n et al.l l2006l. bl ack stars); 13 hr llSevmour et al.ll2008l . orange diamonds); Lockman 
Hole lllbar et al.ll200^ , red circles); and CDF-S l|Padovan^taLll2009l . green upside down triangles). 


with the removal of point sources from the map, extended 
sources were removed also. We also corrected for the re¬ 
duced coverage because of the footprints of observations at 
other wavelengths and for removing radio sources with in¬ 
accurate/contaminated photometry. We find that there is 
a sharp decrease in the overall source counts with decreas¬ 
ing flux density and a flattening below 1 mJy. The AGN- 
powered sources dominate the source counts at bright flux 
densities and their counts decrease sharply at fainter flux 
densities with the form of a power law. The majority of 
other sources emerge only below ~ 1 mJy. 


5 DISCUSSION 

The first thing to note is that many radio sources in the 
ECDF-S contain an AGN. Although this is expected at high 
flux densities, bright radio sources are relatively rare and so 
this result may be surprising for a sample dominated by 
faint sources. The IRAC colour AGN indicator and the cri¬ 
terion for an AGN component being required in the IR SED 
were both conservative. This has resulted in more reliable 


AGN detections compared to a scenario where these detec¬ 
tion thresholds were lower. Furthermore, the effect of the 
FIRRC AGN indicator on the overall detection fraction was 
examined by forcing this indicator to be more conservative 
also. The effect was small. To compare individual indicator 
results to those of previous studies, the fraction of AGN in 
the sub-sample found through the radio spectral index in¬ 
dicator is 22 per cent (13/60 sources with derived spectral 
in dices). This is significa ntly lower than the fraction found 
by ISevmour et all (l2008l) of 69 per cent (31/45). However, 
this is because there was no requirement for the full 3 a con¬ 
fidence range of the spectral indices to lie outside the range 
for the SFGs in the lSevmour et ^ ll2008l) sample. Further¬ 
more, the fraction of sources found to be an AGN through 
their X-ray e mission (52 per cent) is comparable to the frac¬ 
tion found bv IPadovani et ^ (|2009l . hereafter P09) of 58 per 
cent for radio-selected sources, using the same X-ray lumi¬ 
nosity cut. 

We next discuss the IRAC colour distribution of the 
sub-sample. For most of the sources whose radio emission 
is powered by an AGN, it can be seen that they have 
blue IRAC colours which are consistent with old stellar 
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Figure 11. The Euclidean-normalised 1.4 GHz source counts for the radio sources in the ECDF-S. Colours and symbols are as follows: 
total, blue squares; SF-powered, black stars; AGN-powered, red circles; hybrid, orange diamonds; ‘no-z\ green upside down triangles; 
and LERGs, purple triangl es. For clarity, some c ounts have been offset slightly along the abscissa. Also given are the source count models 
from the SKA simulations llWilman et ahllioOSll . 


populations. These objects are likely to be low-luminosity 
low-excitation radio galaxies (LERGs), hosted by ellip¬ 
ticals whose non-s tellar light is weak llLacv et al.l l2004l : 
iGtirkan et akll^Oldh . LERGs do not have a strong ionising 
continuum and therefore they lack the hot dust required 
to pro duce the near - to m id-IR power law SED described 
by the lDonlev et al] (l2012l l criteria. On the other hand, the 
AGN-powered sources that do have a power law SED are the 
more conventional powerful radio-loud AGN such as quasars 
and broad-line radio galaxies (BLRGs) and the radio-quiet 
AGN which have enough radio emission from the AGN to lie 
off the FIRRG. Sources whose radio emission is powered by 
star formation predominantly have slightly redder colours 
and their place in IRAG colour parameter space is_in agree¬ 
ment with that found bv iHatziminaoglou et al.1 (l2009l ). 

The hybrid sources, by definition, are galaxies which 
host an AGN that does not clearly power the radio emis¬ 
sion and therefore they consist of radio-quiet AGN. They 
have the most diverse IRAC colour distribution. The hy¬ 
brids that cover the colour space associated with mid- 
IR-selected AGN are typical QSOs and Seyfert galax¬ 


ies which have a strong ionising continuum, while those 
that have IRAC colours which are similar to SF-powered 
sources are likely to be less powerful radio-quiet AGN. A 
hnal region in the IRAG colour space of hybrid sources 
is one of blue S',/(5.8/rm) / Si. (3.6 pm) and very red 
Si. (8 pm) / S,/(4.5 pm) colours. This region typically con¬ 
sists of low-redshift obscured AGN with strong PAH features 
from a dominant host galaxy dLacv et al.ll200j : iFeltre et al.l 
I 2 OI 3 II . It is reasonable that we hnd that the IRAG colour 
parameter space covered by the AGN-powered and hybrid 
sources extends we ll beyond the region described by the 
iDonlev et al.l (120121 ') criteria. This result is because we have 
identihed many AGN which do not have a dominant AGN 
component in their IR SED, and therefore do not have a 
near- to mid-IR power law. Since the hybrid sources fol¬ 
low the FIRRG, this means their AGN IR emission is sig- 
nihcantly diluted by emission from star-forming processes. 
However, the AGN can still have a n important contribu¬ 
tion to the near- to mid-IR emission (iHatziminaoglou et al.l 
I 2 OO 9 II . particularly for the sou rces that he w i thin t he IRAG 
colour region described by the lDonlev et al] (l2012l f criteria. 





































































Table 4. Euclidean-normalised 1.4 GHz source counts. Columns are as follows: (1) flux density bin range; (2) central flux density of bin in log space; (3) total number of sources in 
bin; (4), (5) (6) and (7) number of AGN-powered, hybrid, SF-powered and ‘no- 2 :’ sources in each bin, respectively; (8) multiplicative correction factor due to the reduced coverage as 
described in Section[4| (9) Corrected effective area, (10) total source counts; and (11), (12), (13) and (14) AGN-powered, hybrid, SF-powered and ‘no- 2 ’ source counts, respectively. 


Si, range 

[mJy] 

(1) 

Central Su 

[mJy] 

(2) 

Total 

(3) 

AGN 

(4) 

No. 

Hybrid 

(5) 

SF 

(6) 

No —2 

(7) 

Correction 

(8) 

AeS 

[deg2] 

(9) 

Total 

[sr-l Jyl'®] 

(10) 

AGN 

[sr-l Jyl-5] 
(11) 

Counts 

Hybrid 

[sr-l Jyl-5] 
(12) 

SF 

[sr-l Jyl'®] 
(13) 

No— 2 : 

[sr-1 Jyl-®] 
(14) 

0.030-0.067 

0.045 

398 

28 

80 

111 

179 

1.30 

0.23 

2.13±0.11 

0.15 ±0.03 

0.43 ±0.05 

0.59 ±0.06 

0.96 ±0.07 

0.067-0.150 

0.100 

139 

32 

28 

32 

47 

1.50 

0.21 

2.60 ±0.22 

0.60 ±0.11 

0.52 ±0.10 

0.60±0.11 

0.88 ±0.13 

0.150-0.333 

0.223 

58 

12 

14 

8 

24 

1.36 

0.24 

3.26 ±0.43 

0.67±0.19 

0.79 ±0.21 

0.45 ±0.16 

1.35 ±0.28 

0.333-0.742 

0.497 

22 

7 

8 

2 

5 

1.45 

0.22 

4.39 ±0.94 

1.40 ±0.53 

1.60 ±0.56 

0.40 ±0.28 

1.00 ±0.45 

0.742-1.653 

1.108 

12 

6 

2 

0 

4 

1.42 

0.23 

7.75 ±2.24 

3.87 ±1.58 

1.29 ±0.91 

0.00 

2.58 ±1.29 

1.653-3.680 

2.466 

6 

3 

0 

0 

3 

1.67 

0.19 

15.14±6.18 

7.57 ±4.37 

0.00 

0.00 

7.57 ±4.37 

3.680-8.194 

5.491 

5 

4 

0 

0 

1 

1.20 

0.27 

30.19 ±13.50 

24.15 ±12.07 

0.00 

0.00 

6.04 ±6.04 

8.194-18.244 

12.227 

1 

0 

0 

0 

1 

4.00 

0.08 

66.86 ±66.86 

0.00 

0.00 

0.00 

66.86 ±66.86 

18.244-40.623 

27.224 

3 

2 

0 

0 

1 

1.67 

0.19 

277.68 ±160.32 

185.12 ±130.91 

0.00 

0.00 

92.56 ±92.56 

40.623-90.450 

60.616 

2 

1 

0 

0 

1 

1.00 

0.32 

369.04 ± 260.95 

184.52 ±184.52 

0.00 

0.00 

184.52 ±184.52 
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Regarding the star-forming total IR luminosities of 
the sub-sample, the greater fractions of NIRGs and LIRGs 
over ULIRGs is simply down to more luminous sources 
being rarer. We compare th e fractions of NIRGs, L IRGs 
and ULIRGs to those for the lSvmeonidis et ahl ll2013ll sam¬ 
ple of SF-powered, IR-luminous galaxies at similar red- 
shifts ( 2 < 2 ). The fraction of NIRGs in our sub-sample of 
43 pe r cent is higher than that for the ISvmeonidis et aP 
(l2013tl sample (20 per cent). We also find lower fractions 
of LIRGs and ULIRGs (53 per cent and 4 per cent, re¬ 
spectively) t han the fractions of 64 per cent and 20 per 
cent for the ISvmeonidis et all ll2013l l sample. This is rea¬ 
sonable because our sub-sample is not IR-selected and some 
sources have been selected because of their radio emission 
from an AGN. Furthermore, SF-powered sources must be 
highly star-forming to be detected in the radio survey, but 
this is not the case for AGN-powered sources. It is be¬ 
cause of this selection effect that higher SFRs are observed 
in the SF-powered sources compared to the AGN-powered 
sources. The low SFRs that we derive for the AGN-powered 
sources are in agreement with those that have been observed 
in other low-redshift AGN hosts ((SFR) = 18—41 M© yr”*^, 
0.4 < 2 < 0.9, ISevmour et ahllioill L 

The fraction of sources powered by star formation that 
host an AGN has been observed to decrease with redshift. 
At z < 0.5, the majority of star-forming sources in the sub¬ 
sample host an AGN while at 2 > 0.5, only a minority of 
sources do. However, this observation is a result of the fact 
that for high-redshift sources, the IRAC photometry were 
not used during the SED fitting because they measure the 
emission at < 2 /rm (rest-frame). Without the IRAC photom¬ 
etry, the AGN component cannot be well constrained and 
thus only low-redshift AGN are identihed by the SED htting 
technique. If the results of the SED fitting are not consid¬ 
ered then the fraction of star-forming sources at 2 ; < 0.5 that 
host an AGN drops to 0.05. 

Regarding the 51 R distribution of the hybrids, there was 
no significant excess in the number of sources with a 51 R 
below the mean 51 R of the SE-powered sources. While on 
an individual basis the radio emission of a hybrid source 
may be AGN-powered, the overall gm distribution indicates 
that the majority of hybrids are SF-powered in the radio. 
Since the hybrid sources are in essence radio-quiet AGN, 
this result supports previous evidence for the radio emission 
from radio-quiet AGN being predominantly a result of star¬ 
form i ng pr ocesses fe.g. iKimball et ^l 201 ll : IPadovani et al.1 
l201ll.l2014lL Focusing on the sub-group of hybrids with an 
AGN identihed through their X-ray emission, Ending no cor¬ 
relation between excess radio and X-ray luminosities is fur¬ 
ther evidence that the AGN in the hybrids is not the source 
of their radio emission. The excess radio luminosity distri¬ 
bution is likely to be a result of residual scatter after the 
subtraction of the emission that is predicted by the FIRRC. 
Studies into the radio and IR emission of radio-quiet Seyfert 
galaxies (e.g. IRov et ^Il998lj have found them to follow the 
FIRRG. Seyfert galaxies constitute some of the X-ray lumi¬ 
nous hybrids and therefore our result is consistent with these 
previous findings also. 

The flux density at which a flattening of the 
Euclidean-normalised source counts is observed is com¬ 
parable to that for the source counts of other sur¬ 
veys, e.g. the IShi XMM-Newton/Chandra Deep Field 


dSevmouretjaL 2008jL Suhaiu/XMM-Newton Deep Field 
(SXD F, Simpson et al.l l2006lj . Lockman Hole (llbar et al.l 
I 2 OO 9 II and a smaller area of the CDF-S covered by P09. The 
total source counts of the EGDF-S are in agreement with the 
counts from ISimpson et al.1 (l2006l l and P09. Below ~0. 2 mjy 
they are slightly lower than those from ISevmour et al.l (1200811 
and llbar et aP (I 2 OO 9 I I. It has been shown that while the 
counts from various radio surveys are consistent at most 
hux densities, below 1 mJy there are some discrepancies 
dSimpson et al.1 12 OO 6 I : ICondon et al.1 12 OI 2 I : iHevwood et al.l 

I 2 OI 3 II . 


Along with calculating the radio source counts for other 
fields, other studies have also separ a ted th e counts by galaxy 
type. For example, ISevmour et al.1 (I 2 OO 8 I I separated the ra¬ 
dio sources in the 13 hr held into SFGs and AGN, although 
they only used radio indicators and so they did not differen¬ 
tiate between SFGs and sources that contained an AGN that 
did not power the radio emission. P09 on the other hand, did 
use non-radio indicators and calculated a radio-quiet compo¬ 
nent of their AGN source counts. Excluding this radio-quiet 
component, P09 found an AGN contribution to the source 
counts of ^10 per ce nt in their faintest bin (43—75/iJy). 
ISevmour et al.l (I 2 OO 8 I I found a slightly higher AGN contri¬ 
bution of ~30 per cent in their faintest bin (38—56^tJy). 
We hnd a contribution of 10—50 per cent from the AGN- 
powered sources in the fainest bin of this study (30—67 /iJy), 
consistent with both studies above. The uncertainty in this 
contribution incorporates the uncertain nature of the ‘no- 2 ’ 
sources in the full-sample. We hnd a contribution of 30—70 
per cent from the SF-powered sources to the counts in our 
faintest bin, in agreement with the contribution from POO 
of ~60 per cent. Regarding the hybrids in this study, we 
hnd they contribute 20—65 per cent in our faintest bin, 
which is consistent with the radio-quiet AGN contribution 
of 15—40 per cent from POO. The contribution we hnd for 
the hybrid sources is also in agreement with the fraction 
of radio-quiet AGN found by iBonzini et al.1 1 I 2 OI 3 II for the 
same radio observations. Increasing the FIRRC AGN in¬ 
dicator threshold to 5 cr ehects the hybrid contribution by 
< 5 per cent. Therefore, increasing this threshold has only a 
small ehect on this result. At higher hux densities, the hybrid 
contribution remains approximately constant up to ~lmjy 
and the hybrid and P09 radio-quiet AGN source counts are 
consistent at 3 ct. Using X-ray data. I^mpson et al.l (l2006lj 
also found that radio-quiet AGN are a signihcant fraction 
of faint (Si, 4 < 300 fJ,3y) radio sources. If we sum the SF- 
powered and hybrid components, the total contribution to 
the counts of 50—90 per cent in the f aintest bin is in agree - 
ment with the SFG contribution from ISevmour et al.l (l2008lj 
of ~70 per cent. Overall, our hndings suggest that galaxies 
which host an AGN but whose radio emission is likely pow¬ 
ered by star formation play a signihcant role at faint hux 
densities. Furthermore, they have a signihcant contribution 
to the battening of the source counts that is observed for 
the extragalactic radio source population. 

Comparing our source counts to the SKA simulations 
as shown in Figure [TT] we hnd our total counts are consis¬ 
tent with those predicted at 1.4 GHz. We hnd that LERGs, 
the AGN-powered objects which do not have a strong ion¬ 
ising c ontinuum and therefore lie outside the lDonlev et aH 
(I 2 OI 2 II region in IRAC colour space, dominate the AGN- 
powered source counts. Using a diherent classification sys- 
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tern, IWilman et alJ (l2008h predict that the less-luminous 
Fanaroff-Riley type I (FR I) class dominates over the more 
powerful type II (FR II) class below SOmJy. A correlation 
between optical emission line and radio luminosities has 
been observed for nearby radio-loud AGN (iKauffmann et ahl 
l2008l l. albeit with a large scatter, indicating that LERGs are 
often FR Is. Our results are therefore in general agreement 
with the radio-loud AGN models from the SKA simulations. 
Furthermore, above O.lmJy the power law slope for the 
LERGs is consistent with the slope for the FR Is. One incon¬ 
sistency between our results and the SKA simulations, how¬ 
ever, is that the hybrid contribution to the counts is approx¬ 
imately half an order of magnitude greater than that pre¬ 
dicted for radio-quiet AGN. With regards to the power law of 
the overall source counts of AGN-powered objects, through 
modelling the evolut ion of the radio luminosity function, 
ISevmour et alj (l2004h showed that this extends to flux den¬ 
sities below ~ 50 /rJy. The SFG component, however, is much 
flatter in the ~50^Jy—ImJy range and the extrapolation 
of the SFG model predicts that at a flux density of ~10 fiJy, 
the SFG source counts will be greater than the AGN counts 
by an order of magnitude. To fainter flux densities, these 
counts will then decrease in a similar manner to the way 
AGN counts de crease from ~1 0 0mJy due to the expansion 
of the Universe. IWilman et al.l (120081 1 and iBethermin et al.l 
(l2012ll both showed that a significant fraction of the SFG 
population at ~ 50 pJy are secularly-evolving galaxies. Note 
that a radio-quiet component was not included in the lat¬ 
ter source count models. Future studies which probe the 
sub-50 fiJy radio source population are required in order 
to test these models at these flux dens ities and to inves¬ 
tigate the nature of even fainter AGN. I Jar vis fc Rawlin'5 
l|2004h modelled the evolution of the radio-quiet component 
of the radio luminosity function and found that radio-quiet 
AGN begin to significantly contribute to the source counts 
at ~ 0.5 mjy. A prediction w hich has been conh rmed by 
later observational studies fe.g. lSimpson et'^l2006ll . Future 
studies modelling the radio counts for the different source 
populations below ~50/rJy will need to present a scenario 
which depicts a strong contribution from radio-quiet AGN. 


6 CONCLUSION 

We have fitted the IR SEDs of 380 radio sources in the 
ECDF-S at 2 <1.4 with models that have ISM and AGN 
components. This has enabled us to separate the IR emission 
from star formation and an AGN in these sources. With the 
additional use of various multi-wavelength AGN indicators, 
our results have revealed that the majority of sources contain 
an AGN. Many of the AGN indicators used in this study 
are conservative and therefore we can be confident that the 
AGN identifications are reliable. We determined that for 25 
per cent of these sources, the radio emission is powered by an 
AGN. Eor 40 per cent, there is no sign of an AGN and their 
radio emission is therefore being powered by star formation. 
The remaining 35 per cent are hybrids in the sense that 
they contain an AGN but their star formation is likely to be 
powering the radio emission. The implication of these results 
is that radio sources which have likely been selected on their 
star formation have a high AGN fraction. 

From the ISM component of the fitted SEDs, some high- 


redshift sources have total IR luminosities that correspond 
to extreme SFRs (> 1000 Mq yr“^). Such high SFRs are 
also observed in Herschel-selected and non-Herschel-selected 
galaxies at higher redshifts than in this study. For the most 
part however, the SFRs of the sub-sample are much lower 
than 1000MQyr“^, with a median SFR of 2OM0yr“'^ and 
a median absolute deviation of 20 Mq yr“^. 

AGN-powered sources dominate the Euclidean- 
normalised source counts at radio flux densities above 
1 mJy. At fainter flux densities, a flattening of the total 
source counts is observed and the contributions from 
SF-powered and hybrid sources become significant. At the 
faintest 1.4 GHz flux density bin in this study of ~50/rJy, 
we find that the hybrid sources contribute 20—65 per cent 
to the total source counts. Therefore, they have a significant 
contribution, along with star-forming sources which do 
not host an AGN, to the flattening of the counts that is 
observed below 1 mJy for the extragalactic radio source 
population. 
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